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Abstract 

Primordial tangled cosmological Magnetic Fields source rotational velocity 
perturbations of the baryon fluid, even in the post-recombination universe. 
These vortical modes inturn leave a characteristic imprint on the temperature 
anisotropy of the Cosmic Microwave Background (CMB), if the CMB pho- 
tons can be re-scatterred after recombination. Observations from WMAP in- 
dicate that the Universe underwent a relatively early re-ionization (z r i ~ 15), 
which does indeed lead to a significant optical depth for re-scattering of CMB 
photons after the re-ionization epoch. We compute the resulting additional 
temperature anisotropies, induced by primordial magnetic fields in the post- 
recombination universe. We show that in models with early re-ionization, a 
nearly scale-invariant spectrum of tangled magnetic fields which redshift to 
a present value of Bq ~ 3 x 10~ 9 Gauss, produce vector mode perturbations 
which in turn induce additional temperature anisotropy of about 0.3 to 0.4 
fiK over very small angular scales, with I upto ~ 10000 or so. 

98.62.En, 98.70.Vc, 98.80.Cq, 95.30.Qd 
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I. INTRODUCTION 



There are two possible processes which might explain the origin of large scale cosmic 
magnetic fields. Both have, however, potential difficulties. One possibility is that some 
high-energy process in the early universe (like inflation or a cosmological phase transition) 
generated primordial magnetic fields which manifest today as galactic or cluster fields [1]. 
The problem here is that this involves speculative physics and there is as yet no compelling 
mechanism to generate fields of the required strength [2-4]. Alternatively, seed magnetic 
fields can get amplified by a large scale dynamo, to produce fields as observed today, coherent 
over galactic or cluster scales [5,6]. There are, however, constraints from helicity conservation 
and/or suppression of lagrangian chaos, due to which the efficacy of this process is unclear 
(see for example [6] and references therein). The effects of weak primordial magnetic fields 
(whose strength today is of order 10~ 9 Gauss), and which are tangled on galactic scales can 
affect galaxy formation [7-11] Hence, it is of considerable interest, to find different ways of 
limiting or detecting such primordial fields (see [2,3,12] for reviews). 

In several earlier papers the consequences of tangled primordial cosmological magnetic 
fields on the observable signatures on the CMBR anisotropy and polarization have been 
investigated in a cosmological scenarios with no re-ionization [13-18]. A detailed numerical 
investigation of CMB signals due to tangled magnetic fields has also been undertaken in 
[19]. The First year results from WMAP satellite observations however indicate that the 
universe could have undergone an early stage of re-ionization [20,21], with a large optical 
depth k ~ 0.17 to the re-scattering of CMB photons. In this paper we follow up our earlier 
work by investigating the consequences of a tangled primordial magnetic field, for CMB 
anisotropics, in a scenario with early re-ionization. In particular we focus on the rotational 
perturbations, that are produced by primordial fields in the post-recombination universe, and 
the additional CMB anisotropy signals that they induce in a universe with early re-ionization. 
Note that compressive velocity modes could also induce CMB anisotropies; however they 
also have non magnetic sources and more importantly suffer larger cancellation effects due 
to the thickness of the last scattering surface around the re-ionization epoch. Hence, our 
focus is only on rotational modes here. 

In section II the general formulation of the problem and the parameters of re-ionization 
relevant for our calculation are derived. Semi-analytic estimates of the additional CMB 
signals are made in section III. In Section IV we present numerical calculations and in 
Section V we discuss our results. 

II. GENERAL FORMULATION 

The equations for the evolution of temperature anisotropy for scalar, vector and tensor 
modes have been derived by Hu and White [22] (hereafter referred to as HW97). We con- 
centrate here on the additional contributions which arise in a re-ionized universe, due to the 
vector modes induced by inhomogeneous magnetic fields. From equation (74) and (56) of 
HW97, the angular power-spectrum of CMB anisotropy corresponding to vector modes is 
given by, 

d-^^fljf-r^w^ap, (2,) 
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Here, V(k, r) is the magnitude of the vorticity, generated by tangled primordial magnetic 
fields, in Fourier space. Also, k is the co-moving wave number, r is conformal time, Tq its 
present value, and ji(z) is the spherical Bessel function of order I. We have ignored a small 
polarization correction to the source term and also a metric perturbation term which are in 
general sub-dominant (cf. [13,16]). The 'visibility function', g(r ,r), is given by, 



where <?(ro, r)dr is the probability that a photon that reaches us at epoch r was last scattered 
between the epochs (r, r + dr). We assume a flat universe throughout, with a total matter 
density parameter Vt m and a non-zero cosmological constant density parameter Vt^ today. 

The exact form of the visibility function is determined by the ionization history of the 
universe. For standard recombination g(ro, r) has only one peak around the epoch of recom- 
bination. However re-ionization can modify this further by making g significant after the 
universe has been re-ionized, the exact modification depending on the complex ionization 
history of the baryons. The modifications of the visibility function will show up in the power 
spectrum of the anisotropy. In this paper we will not be focusing our attention on the physi- 
cal mechanism responsible for re-ionization. We will, however, assume that whatever be the 
source, it results in a sharp transition to a re-ionized situation. We assume this transition 
to be a step function at a redshift z„ corresponding to a conformal time r r {. The value of 
r ri is fixed in this simple model by the optical depth indicated by the WMAP results. The 
visibility function will then have two peaks, one at r rec and the other at r ri . 

The free electron number density in the era after standard recombination is modelled as, 



where ribo = ^b^H^/ (8nGmb) with fit, being the density parameter of the Baryons. (We 
have neglected the small residual electron density after recombination). The optical depth 
is given by 



Here the proper time t and conformal time r are related by using dt = adr and (to) r o) 
correspond to the present epoch. The integration over t can be converted into an integration 
over redshift by the substitution, dt = da/ '(Ha) = —(H(z)(l + z))^ 1 dz. Here H(z) is the 
Hubble expansion rate and is given by 




(2.2) 




(2.3) 




(2.4) 



H{z) = H [Q A + Q m {l + z) 3 ] 1 / 2 . 



(2.5) 



The optical depth up to re-ionization can now be expressed as, 




(2.6) 



where y is defined as y = (1 + z) 3 (Q m /Q\). On integration, this gives, 
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(2.7) 



Adopting, K ri = 0.17, h = 0.71, fi b = 0.044, Vt m = 0.27 and tt v = 0.73, we get, 
z ri = 14.57. (This value is close to the one obtained from equation 24.79 given in Peebles 
[23]. We have neglected a small correction due to helium fraction. This, however, does not 
significantly affect the results.) The conformal time of re- ionization is given by, 



(yja ri + a eq - y/a^j (2.8) 



where, a ri x = 1 + z ri and a 1 = 1 + z eq specifies the scale factor at matter-radiation equality. 
The current epoch is given by 

t = [1 - 0.0841 ln(n m )} (yl + a eq - v ^) (2.9) 

With 2CHQ 1 = 6000 h _1 Mpc, z eq = 3233 and z ri = 14.57, we get r = 12595fc _1 Mpc, 
r r j = 2741/i _1 Mpc and ro — r r j = 9854/i _1 Mpc. Using the functional form of the number 
density of the electrons as given in equation (2.3) and the r ri determined above, we can 
calculate the form of the visibility function. We will use this form for the visibility function 
in the numerical computation in section IV. However, to begin with, we approximate the 
visibility function ^2(^0, r ) after recombination, to be a truncated exponential, and estimate 
the CMB anisotropy in a semi-analytic manner. 
Specifically, we adopt, 

92(r ,r) = ^l e - T -^Q{T-r ri ) (2.10) 
a 

Here the Heavyside 6(x) function, is zero for x < and 1 for x > 0. It takes account of 
the fact that before re-ionization, n e is negligible. Further, N2 is a normalization consant 
and a gives the spread of the exponential. By appropriately choosing a, we can set the 
width of the reionized last scattering surface. Also note that g(r , r) has the interpretation 
of probability; so its integral over r from r = to r = ro should be normalized to unity. 
This determines the normalization factor N 2 . For a sufficiently early epoch of re-ionization, 
we generally have (r — r ri )/a ^> 1. In this case, the condition that the integral of g(r ,T) 
over r should be unity implies N 2 + e~ Kri = 1, or N 2 = 1 — cxp —(n ri ). So N 2 measures the 
probability of at least one scattering between ro and r ri , due to the re-ionization. For small 
K r ; < 1, we have N 2 ~ n ri 

The constant a is determined by the (conformal) time, say r m , after which g 2 drops to 
1/e times its peak value (at r ri ). Thus if 



02(Tp,T m ) = 1 

5-2(7-0, T ri ) e 



(2.11) 



then, a = r m — r ri . To determine r m , we use the exact form of g 2 and calculate the epoch 
when [g(r ,T ri )/ g(r ,T m )} = e. Using equation (2.3) and the expression for the visibility 
function in equation (2.2), we get, 
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9(TQ,T r i) 

9(ro, T m ) 



n P (Tri a(T r . 



n e (r m )a(T m ) 



cxp 



n e (T")a T a(T")dT'' 



(2.12) 



In the interval between r ri and r m , the universe is in general matter dominated. Hence 
a(r) oc r 2 . Also for r > r ri , we have, n e oc a~ 3 . Hence, we can simplify the above equation 
to 



cxp 



—C(J T - 



f 



8nGrrib 

Substituting for the scale factor as a = t 2 /tq, we get, 

, 'TrrA _ OmO Hq Of CT r j ( 1 + Z ri ) 2 



(2.13) 



87iGmt 



(2.14) 



This gives r m = 3519/i _1 Mpc for r ri = 2741/i~ 1 Mpc giving the value of a as 778/i~ 1 Mpc. 



III. RESULTS OF SEMI- ANALYTIC APPROXIMATION 



Consider first the limiting case when ji is more sharply peaked than g^. This implies 
the limit I ^> (r — r ri )/a. On using the values of r , r ri and a determined above this limit 
translates to I 3> 12. We will also assume that the source term V varies slower than either 
the variation of g<i or that of ji. These approximations will help in understanding the effect 
of the model with re-ionization in a semi-analytic manner. Results of numerical calculations 
which do not make these approximations are given in the section IV. 

Since we wish to focus on the effects of re-ionization we take the visibility function to be 
g 2 in the integral appearing in equation (2.1), which then becomes, 

Jo dr 92(r ,r)V(k,r) ] ^— ) . (3.1) 

The function, ji peaks at / ~ k(r — r) or r ~ r — l/k. So in the above integral we evaluate 
V(k, t) and ^(^Oi r) at r = r — l/k and replace k(r — T ri ) by / and move them out of the 
integral. We also use the identity, 

jf aw «fc = Vs- < 3 - 2 > 

Substituting the final result in the expression for Ci in equation (2.1), we get, 

Q = jT dk k 2 ( | V(k, r - l/k) | 2 > e -^-w-i/*)/a e(To _ Tn _ l/k) (3. 3) 

t\t2 r ( To - T A 

= tAj dx-( | V(k,r - l/k) | 2 }e-^Q(x) (3.4) 

2a z l Jo I 

where, we have substituted x — {tq — r r i — |). The exponential term peaks at x — 0. So we 
will evaluate the rest of the integrand at x — > or k — > Z/(r — r ri ) and move it out of the 
integrand. The integration of the exponential factor then gives, 
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c. 



N 2 2n 2 /Tq-t, 



AaP V I 



k 3 ( | V(k, t — l/k) 
2~7 2 



(3.5) 

k=l/(T -T ri ) 



Kl + l) -Q = ^ f ^ ) A 2 v(^ (3.6) 



2tt 4a I i 

for I >> 1 and r — r ri ^> a. Here we have defined the power spectrum associated with 
rotational velocity perturbations, 

K(k , T) = iumim . ( , 7) 

We assume that the magnetic field which induces vortical perturbations is initially a 
Gaussian random field. On large enough scales, the induced velocity is generally so small 
that it does not lead to any appreciable distortion of the initial field [10,24,25]. So, the 
magnetic field simply redshifts away as B(x, t) = B (x)/a 2 . The Lorentz force associated 
with the tangled field is then F L = (V x B ) x B /(47ra 5 ), which pushes the fluid and 
creates rotational velocity perturbations. These can be estimated as in [8] by using the 
Euler equation for the baryons and we give a detailed derivation in Appendix A. 

Further, in order to compute the ensemble average Ay and hence the Qs, we need the 
magnetic spectrum M(k). This is defined using < &j(k)6j(q) >= 8^ tCl Pij(\s)M(k), where 5k, q 
is the Kronecker delta which is non-zero only for k = q. Here b(k) is the Fourier transform 
of B , the present day value of the tangled magnetic field. Also P»j(k) = ((% — kikj/k 2 ) is 
the projection operator which ensures that magnetic field has zero divergence. This gives 
< Bg >= 2 f(dk/k)Al(k), where Ajj (&) = k 3 M(k)/(2n 2 ) is the power per logarithmic inter- 
val in k space residing in magnetic tangles, and we replace the summation over k space by 
an integration. It is convenient to define a dimensionless spectrum, m(k) = A 2 (/c)/(_Bq/2), 
where B is a fiducial constant magnetic field. The baryonic Alfven velocity, Va, for this 
fiducial field is, 

^ = Ta B \ lf 2 ~ 15 x W 5B -9, (3-8) 

where pbo is the present day Baryon density, and _B_g = (-Bo/10~ 9 Gauss). 

We will also consider as in [15], power-law magnetic spectra, M(k) = Ak n . We will 
cut-off this spectra at the scale where the perturbations are no longer linear, say at k — /cjv- 
We expect k N to be of order galactic scales, or k N ~ (fhMpc)^ 1 with / ~ 1, for the range 
of redsifts which make a non-zero contribution to the visibility function and for the field 
strengths that we consider [10,25]. We fix A by demanding that the field smoothed over a 
scale, kN, (using a sharp /c-space filter) is B , giving a dimensionless spectrum for n > — 3 
of 

m(fc) = {n + 3){k/k N ) 3+n . (3.9) 
Assuming such a spectrum we have from Appendix A, 



k 3 ( | V(k,r) | 2 ) /2c\ 2 fl + z 



2tt 2 \h J \n 



kV 



V8 



I(k) 



(3.10) 
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Here I(k) is a mode coupling integral which has been worked out in detail in [18,13]. In 
particular for the case n < — 3/2 which we consider below it is given by, 

J 2 (A;) = |(n + 3)(A )6 +- (3.11) 
3 Kn 

Using equation (3.11) in equation (3.10) and substituting the resulting expression for velocity 
power-spectrum in equation (3.6) we get, 

'(' + l).G = ^f#V VX 5(„ + 3)^f±) ,+2 " ,3,2) 

N J u— I 



2tt 4 \HJ m m 3 a \k 



We define AT = T [l(l + l)Ci/2ir} 1 / 2 as a measure of the temperature anisotropy, where 
T = 2.73 K is the present CMB temperature. Taking the square root of the above expression 
and using the values of Va, fcjv, Q m and z„ mentioned earlier we get, 

AT = 0.52 N 2 B 2 _ g (n + 3) 1/2 f~ {n+3) (-J—) 3 - 5+n ^ K ( 3 _ 13 ) 

9854 

(Note that k N (r — r ri ) = /9854). For a field of 3 nano-gauss (B_ g = 3), a nearly scale 
invariant spectrum for the magnetic field (n = —2.9) and for an optical depth of ~ 0.17 or 
N 2 ~ .17, we get, 

/ / \°- 6 

AT = 0.25 (3.14) 

\9854 7 v ; 

A similar analysis in the opposite limit (g 2 more sharply peaked than ji an d hence 
I < (r - T ri )/a ~ 12) gives 

AT = 5.3x A^yn + 3) 1 / 2 ^^ pK (3.15) 

For a field of 3 nano-gauss (B_ 9 = 3), a nearly scale invariant spectrum for the magnetic 
field (n = —2.9) and for an optical depth of 0.17 or N 2 ~ .17, we now get, 

AT = 2 - 66 {<mT » K (316) 

For example at I ~ 10 one predicts 1.3 x 10~ 3 /iK, which is very small compared to other 
signals expected at these low / values. 



IV. NUMERICAL RESULTS 

While the calculations in the semi-analytic approximation in the last section provide us 
with rough estimates, to make concrete predictions we need to numerically compute the 
temperature anisotropy. In this section we give the results of the numerical evaluation. 
Specifically, we numerically evaluate the integral in equation (2.1). As discussed earlier, 
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the visibility function appearing in this integral has a dominant contribution at the epoch 
of standard recombination and another at the epoch of re-ionization. The contribution at 
the epoch of re-ionization we have denoted by g 2 and we are interested in the additional 
contribution to C\ resulting from g 2 . Unlike the semi-analytic case, (where we approximated 
g 2 as an exponential decay for epochs later than the re-ionization epoch) we use the form 
given in equation (2.2) with the number the number density of the electrons as given in 
equation (2.3). The resulting expression for g 2 is, 



92(t ,t) 



3Q m0 H$a T Tq 
87rGm 6 r 4 



6(t - T ri ) exp 



P dr' (to/t') 



(4.1) 



Further, as the universe is believed to have been matter dominated after r = r ri , the redshift 
and the conformal time are related by, 1 + z— (t /t) 2 . We can neglect the accelerated phase 
of the universe, as this phase is believed to have set in at a low redshift (z ~ 1) by which 
time, the visibility function would have decayed sufficiently Hence, this will not introduce 
any significant error in our computation. With these simplifying assumptions, we have 
computed AT by evaluating the r and k integrals numerically in Eq (2.1). While evaluating 
this we have retained the analytical expression for I{k) given in equation (3.11). For B ~ 3 
nano Gauss and a nearly scale invariant spectrum (n = —2.9) we find AT = 0.33/iA" for 
I = 10000. For higher n the fluctuations are larger. The results of the numerical calculation 
are shown in Figure 1, for the magnetic spectral index n = —2.9, —2.8 and —2.7. We 
can compare this with the semi-analytic results from equation (3.14). For n = —2.9, our 
semi-analytic calculations give AT = 0.25 for / = 10000. So we see that although the semi- 
analytic calculation marginally underestimates the temperature anisotropy, we do get the 
correct order. This is the observed trend for all n. For n = —2.8 and n = —2.7, AT turns 
out to be 0.47 [i K and 0.58 n K respectively for / ~ 10000. Indeed both the amplitude of 
AT and its / dependence, computed using the semi-analytic calculation, agrees reasonably 
well with the more exact numerical integration. 



V. DISCUSSION AND CONCLUSIONS 

We have investigated the additional CMB temperature anisotropies that are generated 
by tangled cosmological magnetic fields in a Universe that underwent a relatively early re- 
ionization. We are motivated by WMAP results, which indicate that the universe could have 
been re-ionized as early as z ~ 15. We have focused for the present on rotational velocity 
perturbations, which can be substained only by cosmological magnetic fields. These modes 
also suffer a milder damping due to the finite thickness of the re-ionized last scattering 
surface, than compressional modes. Our results supplement earlier work obtained in the 
context of a Universe that does not undergo re-ionization. 

We find that a nearly scale-invariant spectrum of tangled magnetic fields (with n = —2.9 
to n — —2.7), which redshift to a present day value of about 3 nano Gauss can produce 
anisotropies at the level of about 0.3/xK to 0.5yuK for I ~ 10000. Even larger signals would 
obtain if we were to consider models with larger n. We have simply stopped at this large / 
because we cannot use linear theory at present to calculate the expected signals for I larger 
than /^(tq— r ri ) ~ /9854. It is also interesting to note that, at these large /, the above signals 
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are comparable to primary signals due to tangled fields, arising from the usual recombination 
epoch [26]. CMBR anisotropy experiments which probe such very small angular scales can 
thus be useful to measure or impose bounds on the magnitude as well as the spectral index 
of these magnetic fields. We have also considered here only the simplest first order effects 
due to re-ionization. It would be of interest also to estimate the effects of an inhomogeneous 
re-ionization on the CMB anisotropy from vector modes. 



APPENDIX A 

We need to calculate the rotational velocity field of the baryons in the post-recombination 
era, in particular around the epoch of re-ionization. We shall follow the formalism developed 
by Wasserman [8] . Here we summarize briefly the essential features of the derivation based 
on Wasserman's paper. 

The Fourier space, linearized Euler equation for the rotational component of the baryon 
velocity in the post recombination era is given by, 

dt a 1 lnp b {t)a* 1 } 

Here Vi(k,t) is the Fourier component of the rotational velocity perturbation, Fj is the 
Fourier component of the vector [(V x Bo) x Bo]j and Pij = 5ij — kikj/k 2 as usual projects 
out its rotational component. Note that V(k,t) is the magnitude of the vector Vi(k,t), that 
is | V" | 2 = ViV*. We have also assumed here, (as mentioned in the text) that on large enough 
scales, larger than A;^ 1 ~ Mpc, the velocities are so small that the magnetic field does not 
get significantly distorted, but simply redshifts away as B(x, t) = B /a 2 . 

The above equation can be solved in a straightforward fashion to obtain the rotational 
component of the velocity field. Note that in the post-recombination era, we can neglect 
the radiation density compared to the matter density. Further, the dark energy component 
dominated over matter at late times, i.e., at redshift less than unity, whereas the results from 
WMAP can be interpreted to mean that the Universe underwent a re-ionization phase at an 
earlier time. Thus from the epoch of recombination (t rec ), till the epoch of re-ionization (t r i), 
pressureless matter was the dominant component of the Universe. Hence, for t rec < t < t r i 
we can take a(t) = a rec (t/t rec ) 2 ^ 3 = (1 + z)^ 1 . The solution for Eq. Al is then given by 



^t rec (Pij Fj 



A-np m a 



2 

rec 



t \ 1 / 3 ft \ 2 / 3 * 
"rec \ I v rec \ 



(A2) 



Here pbo is the present day Baryon density, and we have assumed that the rotational velocity 
was negligible at recombination. (Any small rotational velocity at recombination contributes 
to a faster decaying term in the above solution than the term we will retain below). Since 
a(t) oc t 2 / 3 , H rec = 2/(3t rec ). Also we have from Einstein equation H 2 ec = 8irGp rec /3 = 
H 2 tt m0 (1 + z rec f. Thus 

2 (1 + z )~ 3/2 

tree — ott /o ^ > 

Substituting this expression for t rec in Eq. A2, and noting that a(t) = 1/(1 + z), we have, 
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Vi = 



2 1 PjjFj 
Ho V^mo 4vrp fe0 



x(l + z) 



1/2 



(A4) 



For a power spectrum of B as given in the text, a tedious but straightforward computation 
gives the power spectrum of the rotational component of the Lorentz force. We have 



e < | PvFjPaF? | ) 



2tt 2 



4np b 



P\/ 4 



(A5) 



where I(k) is a mode coupling integral whose explicit form is given in [10,15]. The magnetic 
field enters through the baryon Alfven velocity Va- The power-spectrum of the V (or the 
magnitude of Vi) required in the text is then given by 



2tt 2 



(\V(k,z) \ 2 )=( W - 



1 + z 



V8 



I(k) 



(A6) 
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FIGURES 

FIG. 1. The figure shows the temperature anisotropy due to vector type modes induced by 
tangled cosmological magnetic fields that are of a strength of 3 nano Gauss today. The magnetic 
power spectrum is assumed to be a power-law characterized by an index n. The figure shows the 
results for n = —2.9 (solid line), —2.8 (dotted line) and —2.7 (dashed line). 
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